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We study the magnetic field dependence of the dielectric 
response of large cylindrical molecules such as nanotubes. 
When a field-induced level crossing takes place, an applied 
electric field causes a linear instead of the usual quadratic 
Stark effect. This results in a large dielectric response. Ex- 
plicit calculations are performed for doped nanotubes and a 
rich structure in the real part of the low-frequency dielectric 
function e (H) is found when a magnetic field is applied along 
the cylinder axis. It is suggested that studies of e'{H,T) can 
serve as a spectroscopic tool for the investigation of large ring- 
shaped or cylindrical molecules. 



PACS. 77.22.-d Dielectric properties of solids and liquids 
78.40.Ri Fullerenes and related materials 
75.20.-g Diamagnetism and paramagnetism 
73.23.-b Mesoscopic systems 



Dedicated to K. Biedenkopf on the occasion of his 70th 
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I. INTRODUCTION 

During the last years considerable progress has been made 
in precise measurements of the real part of the low-frequency 
dielectric function ^'(i^). In particular those measurements 
could be extended to ultra- low temperatures, i.e., down to a 
few mK. Ratios of 5e' /e' up to 10"^ were achieved in that tem- 
perature regime. A rather spectacular success associated with 
that progress was the observation of a strong magnetic-field 
dependence of the polarizability of multicomponent glasses 
in the mK regime This development suggests reconsid- 

eration of magnetic field effects on ring molecules or related 
structures like nanotubes. In both cases an applied magnetic 
field induces a diamagnetic ring current. Due to this cur- 
rent the energy of the ground state increases quadratically 
with the applied magnetic field. This continues until one of 
the excited states which lowers its energy in a field crosses 
the ground-state and becomes the new ground state. At the 
crossing point an applied electric field causes a linear Stark ef- 
fect, instead of the usual quadratic one and hence a divergent 
electric polarizability. The physical origin of the crossover is 
easily understood. An excited state carrying a ring current 
in the absence of a magnetic field becomes a state without 
a ring current, when a sufficiently high magnetic field is ap- 
plied, because the induced current may cancel the original 
one. When this is the case the energy of that state equals 
the one of the ground state in the absence of a field. This 
simple argument shows that the ground-state energy is a pe- 
riodic function of an applied magnetic field. The periodicity is 
given by the fiux enclosed by the ring current. When one flux 
quantum (j>o = hc/e is penetrating the ring the ground-state 
energy has returned to its original value. It is well known 
that a huge field of order 10^ T is needed for a fiux quantum 
(po to penetrate a benzene molecule consisting of a ring of six 
carbon atoms. Because of that large value of the magnetic 
field, effects of it on ring molecules have obtained only little 
attention in the past The purpose of the present paper 

is to point out that the situation has changed considerably. 
Not only has there been experimental progress in performing 
high precision measurements of e' at low temperatures, but 
also the synthesis of organic ring structures has made impres- 
sive advances. For example, nanotubes of large circumference 
have been produced, lowering the field required for the en- 
closure of a fiux unit. In this paper we want to demonstrate 
that a measurement of e'{H,T) should provide important in- 
formation on the electronic excitations of ring molecules, in 
particular on level crossings. As a first step, we calculate here 
the dielectric response of a molecule consisting of a square lat- 
tice, e.g., of carbon sites bent into a cylindrical form and of a 
nanotube, i.e., a bent honeycomb lattice. Various extensions 
of the work presented here will follow later. 



1 



The magnetic and electric field are assumed to be directed 
along the cylindrical axis. The electron interactions are as- 
sumed to be included in effective one-electron parameters like 
in an extended Hckel theory or in the quasiparticle theory of 
Landau. In a subsequent investigation we shall include the 
electronic interactions more explicitely than done here. This 
may have profound effects on the results. A rich structure 
in e'(H,T) is obtained which should be experimentally ob- 
servable. It is closely related to the low-energy excitations 
of the systems in an applied magnetic field and in particular 
to level crossings as the field changes. Although our find- 
ings are limited here to the cylindrical structures described 
above, they suggest detailed experimental studies of e'{H,T) 
for ring molecules. This seems to us a largely untouched field 
of research. 

In order to explain the main features of e'{H, T) we consider 
first a single ring of N sites in a magnetic field along the ring 
axis. The Hamiltonian is 



H ■ 



-Fh.c.) 



(1) 



where a^^,ana are electron creation and annihilation opera- 
tors and (j> is the magnetic flux through the ring in units of 
the flux quantum 00 . The resulting energy eigenvalues are 



2-K , 
"/V 



:0,±1,±2,.. 



(2) 



The ground-state energy Eg is periodic in the flux, i.e., 
£'g(0-|- 1) = Eg{<p). More explicitely we write for Eg 



-sin(^g) sin(^0) 



(3) 



For a closed-shell system, i.e., for an electron number iVg = 
4n + 2, where n is an integer we find that 



E. 27r 
sm—q = 



(4) 



In that case the field-dependent contribution to the ground- 
state energy is 

SEg{^) = Y^2t cos (^g) (cos (^0^ - l) 

occ 

= -£;,(O)(l-cos(^0)) . (5) 



For large N we may expand this expression and obtain 

2n^ 



5Eg{4>) = -EgiO) 



JV2 



> 



(6) 



for < i . When <f) = \ the ground state is twofold degenerate 
because of a level crossing at that point. For > i, the 
expression (3) is replaced by 



5Eg(<f>) = -Eg{Q) ^ (1 - c^f 



(7) 



The behaviour of SEg{(f>) is schematically shown in Fig. 1. 
The contribution of 5Eg{(j)) to Eg{(j)) is very small for large 



values of N and hardly detectable. This does not hold true 
though for other quantities. For example, when in addition 
an electric field is applied perpendicularly to the ring axis its 
effect on the ground-state energy is strongly dependent on cf). 
For (/) = i we are dealing with a linear Stark effect instead off 
a quadratic one when <^ ^ i . Therefore the dielectric constant 
has a singularity at that particular value of (p. This simple 
example sheds light onto the physical reason why the dielec- 
tric function can be so sensitive to an applied magnetic field. 
The same feature is found for cylindrical molecules which are 
subject of this paper. 



II. MAGNETIC FIELD DEPENDENCE OF THE 
FREE ENERGY 

In order to demonstrate the influence of an applied mag- 
netic fleld H on the free energy we consider two different 
systems. One is a model square lattice rolled into the form of 
a cylinder. The other one is a nanotube which consists of a 
honeycomb lattice rolled into a cylinder in the same way. 

We start with the square lattice forming a cylinder. It 
consists of N atoms along the perimeter and of N atoms along 
the cylindrical axis z. The eigenvalues depend on the flux <f) 
though the cylinder and are of the form 



27r 



+2n 



N 

2 me«CT 

m 



(p + <^) 



teos 



2n 



M + 1 
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The first term corresponds to a transfer integral of size — i 
along the perimeter and the second to one of magnitude — | 
along the z axis. The parameters p and q take the integer 
values p = 1,...,N and q — 1,...,M, respectively. The last 
term is the Zeeman contribution which is expressed here in 
terms of the flux 0. Since the latter is in units of the flux 
quantum 0o, the ratio of the effective mass meg divided by the 
electron mass rn enters here, with the former referring to an 
electronic motion perpendicular to the z axis. Furthermore, 
cr = ±1. The free energy of the system is of the usual form 



I3F=-Y^ ln[l + exp(-/3(ep 



(9) 



where P = (fcsT)"^ and /i is the chemical potential. It is 
determined by expressing the number of electrons A^e in terms 
of it, i.e.. 



(10) 



In practice we calculate by first choosing an approximate 



value 110 and calculating the corresponding value Ai^'. The 
correction 5fj,o to /Uo can then be determined from 



1 / A*'"'-' — AT 



(11) 



This expression is more convenient for numerical calcula- 
tions than its linearized version in (Ai^' - Ae)/a(r) where 
the function a{T) is given by 
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a{T) = i cosfe" 



0f 

2 l^p<Jo- 



Ato) 



(12) 



One can use the corrected potential fii = fiQ + Sjio in or- 
der to calculate the next correction Sfii. We obtain with 
5^1 the chemical potential already with an accuracy of order 
N~'^M~^, which perfectly serves our purposes. 

The same procedure can be applied to carbon nanotubes. 
In that case the unit cell contains four carbon atoms. Hence 
the excitation energies form four bands, i.e., 



,(<^) = ± l + ?ip±(l+Up«,)! 
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with 



Up = 2 ^ 



1 + cos 



2. 

N ^ 



w„ = 2 ( 1 + cos 
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and p = 1, , A''; q — 1, . 



.M 



and Eka denotes the excitation energies of the chain. In order 
to compute D from ( [l6| ) we have to diagonalize Hid in order 
to find the eigenenergies and eigenfunctions of that Hamil- 
tonian. This can be done if not more than 1000 atoms are 
involved. Instead of doing that we want to use here a simpler, 
more effective quasiclassical scheme. In the quasiclassical ap- 
proximation the excitation spectrum is of the form 



= — 2cos 
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(18) 



with p = 1, ,M and m — 1, ,A4. The corresponding 

expression for the induced dipole moment is 



D 



2eao 
M 



E 



(m-(M-H)/2) 

g;3(eprn-A') _|_ \ 



(19) 



We have calculated D for a chain of M = 201 atoms by 
using ( ]l^ ) and alternatively (p^. The results are compared 
in Fig. 2 for different densities and temperatures. The devia- 
tions caused by the semiclassical approximation are less than 
1% or -jj in all cases. This justifies the use of a quasiclassical 
approximation when we calculate the dielectric response of 
cylindrical molecules such as nanotubes in an applied mag- 
netic field. 



III. INDUCED DIPOLE MOMENT 



When an electric field is applied along the z axis the ex- 
citation spectrum of the system can no longer be calculated 
exactly. Instead, approximations have to be made. We cannot 
apply linear response theory because of a linear Stark effect 
at level crossings. Since in practise the applied electric field is 
very small, the density of electrons changes only slightly along 
the z axis. This enables us to determine the induced density 
changes by using a quasiclassical approximation. Within that 
scheme the excitation energies depend not only on p, q and a 
but on the coordinate z as well. We illustrate the approxima- 
tion by considering a chain of Af atoms as a simple example. 
The Hamiltonians is of the form 



Hid = - ^ (aJ,+i)„an<T + h.c.) 

n — l,a 

+ eaoFo ^ a^^a„a (n — ^ ^ 



(15) 



Here ao is the lattice constant and Fo is an applied electric 
field along the chain direction. 

Exact calculations of the induced dipole moment D require 
an evaluation of the expression 



D = 



Sp d e-'^'^id-M) 
Sp e-''(^id-M) 



Here d is dipole operator 

M 



M + 1 



(16) 
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IV. RESULTS AND DISCUSSIONS 

In the following we want to present results for the dielectric 
response of the two types of cylindrical molecules described 
above, i.e., for a square lattice rolled into a cylinder and for 
nanotubes. The induced dipole moment is calculated in close 
analogy to the one of a ring, although here the electric field 
-Fo is directed along the cylindrical axis. We start with the 
square-lattice case. In analogy to (^^ the induced dipole is 
calculated from 



(m- (M + l)/2) 

g/3[ep5CT{"i,0)-M(0)] _j_ \ 



(20) 



where 



r(m, (f)) 



,{<f>) + eaoFo (m - (Af + l)/2) (21) 



and €pqa{4>) is given by (|8|). Results for the magnetic-field 
dependent part D(<^) — D(0) are shown in Fig. 3 for a cylinder 
with 100 atoms along the circumference and 1000 atoms along 
the axis, i.e., = 100 and M — 1000, respectively. Note that 
{D{(t>) - D(0))/D(0) = (e'(0) - e'(0))/e'(0) where e' is the 
real part of the dielectric response in the low-frequency limit. 
The temperature, or more precisely ksT is 10~* in units of 
the hopping matrix element. The chosen density corresponds 
to 0,74 electrons per site. One notices a rich structure as a 
function of the applied magnetic field. 

A cylinder formed from a square lattice is a hypothetical 
case. But if one assumes a lattice constant ao — 1.40A as in 
the case of an aromatic carbon ring, the field required for the 
enclosure of a flux unit is of order 260r. It is derived from 
the following relation between the flux <j) (in units of (po) and 
the applied magnetic field 
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N^al eH 



(22) 



The structure in e'(-ff) obtained within that range of fields 
reflects properties of excited states, in particular crossings of 
energy levels. 

For nanotubes the calculations are done quite similarly, but 
here we have to take a sum over all four energy bands. The 
computational results are shown in Fig. 4 for a density of n = 
0.89 TT-electrons per site. One notices that the rich structure 
in 5e'{H)/e'{H) in the regime < H < AOT is of order unity 
and therefore should be easily detectable. Results for other 
densities look similar, except for n = 1 which is special. The 
reason is that a honeycomb or graphite lattice has for n = 
1 a Fermi surface consisting of a point. Therefore in a a 
flnite system the level spacing close to the Fermi energy is 
particularly large. This leads to small changes in 5e'(_ff)/e'(0) 
only. The situation changes at high magnetic fields. Due 
to the Zeeman term in the Hamiltonian the spin dependent 
densities ria differ more and more from each other, i.e., n — 
rij + with n^{H) ^ ni{H) and the Fermi surface moves 
away from the special point at half filling. This brings us 
back to the doped case and we obtain again a rich structure 
in 5e[H)/e'{0) like in Fig. 4. 
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V. CONCLUSIONS 



The above calculations show that large molecules of cylin- 
drical or circular shape should show detectable magnetic fleld 
effects due to the Bohm-Aharonov effect. They lead to a 
strong variation of the dielectric function in the low-frequency 
limit as function of the applied magnetic field. Those varia- 
tions are predominantly caused by doubly degenerate ground 
states resulting from level crossings in the applied field. At 
a crossing point an applied electric field causes a linear Stark 
effect instead of a quadratic one when the ground state is 
nondegenerate. The present investigation requires a number 
of extensions which will be the subject of separate investi- 
gations. One concerns the dependence of 5e{H) on the di- 
rections of the applied magnetic and electric fields. Impor- 
tant is also a proper inclusion of electron correlations. As 
pointed out before, the present calculations have been done 
within the one-electron approximation. But correlations, in 
particular when they are strong will clearly result in impor- 
tant modifications of the dielectric response Finally, 
we also have to generalize the above theory to the case of 
mutually interacting molecules. This may become an impor- 
tant issue when platelike molecules are forming stacks and a 
magnetic fleld is applied along the direction of the stack. Al- 
though the work presented here needs extensions of the form 
just described it is fair to state that the results presented 
here justify efforts towards a systematic investigation of the 
magnetic-field dependent dielectric response of ring- or cylin- 
der shaped molecules. We feel that in the future they may 
develop into a spectroscopic tool for studying low-energy ex- 
citations of such systems. 
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FIGURE CAPTIONS 

Fig. 1: Schematic plot of (5i?g{(?!>) </> (thick solid line). The 
level crossings are due to parabolas describing differ- 
ent excited states ei{4>) shifted by flux units. When 
levels cross, an applied electric field causes a linear 
Stark effect and hence induces a divergent dielectric 
response. 

Fig. 2: Induced dipole moment as a function of electrons per 
site for a chain of M = 201 atoms calculated with the 
exact quantum-mechanical expression ( |l6[ ) (dashed 
lines) and when a semiclassical approximation (19[) 
(solid lines) is made, (a) and (b) correspond to tem- 
peratures kgT = 0.01 and 0.05, respectively (in units 
of the transfer integral). 

Fig. 3: Dielectric response [e'{H) — e'(0)]/e'(0) for a model 
square-lattice system with N = 200, M = 1000 in an 
axial magnetic field. The temperature is kgT = 10~^ 
(in units of the transfer integral), and ao = lAA. The 
density is 0.74 electrons per site. 

Fig. 4: Dielectric response [e'{H) — e'(0)]/e'(0) for a nan- 
otube with N = 100, M = 1000 in an axial magnetic 
field. The temperature is fesT = 10~^ in units of 3 
eV and the density is 0.89 electrons per site. 
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